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ABSTRACT 

This paper aims to investigate the influence of inhomogeneous radiation in void for¬ 
mation. Since the process of void formation is non-linear, a fully relativistic model, 
which simulates the evolution of voids from the moment of the last scattering until the 
present instant, is presented. It is found that in order to obtain a model of a void which 
evolves from Sq = —10 -5 at the last scattering moment to the present day S = —0.94, 
the existence of radiation must be taken into account. The ratio of radiation energy 
density to matter energy density in CDM models at the moment of last scattering 
is 1/5. This paper proves that such a value of radiation energy density cannot be 
neglected and influences the first stages of void evolution. Namely, it is important to 
the process of structure formation and hence significantly influences the dynamics of 
the Universe in the first millions of years after the last scattering. 

From the fact that the evolution of voids proceeds differently in various cosmo¬ 
logical background models we use the process of void formation to put some limits on 
values of cosmological parameters. We find that the model with fl ma t ~ 0.3 fits the 
observational data best. 

Key words: Cosmology: cosmic microwave background; cosmological parameters; 
theory; early Universe; large-scale structure of Universe. 


1 INTRODUCTION 

At the end of the 1970s several galaxy redshift surveys 
started to operate. As a result the spatial distribution of 
a large sample of galaxies had been measured. It turned 
out that galaxies are distributed inhomogeneously and form 
structures like voids or clusters. 

The most probable explanation of existence of such 
structures is that they evolved from small initial fluctua¬ 
tions. These fluctuations can be traced in the observations 
of the Cosmic Microwave Background (CMB). 

While gravity is an attractive force, it is relatively easy 
to reproduce high density regions, for instance by setting 
the initial conditions so that collapse or shell crossings oc¬ 
cur. This cannot be done in case of low density regions, 
such as voids, inside which the estimated value of the den¬ 
sity contrast is less than 5 = —0.94 (which is less than 6% 
of the mean background density) (Hoyle & Vogeley 2004). 
The description of evolution of voids from initial fluctuations 
consistent with observational constraints is very difficult. So 
far none of the attempts to solve this problem succeeded. 
Even N-body simulations predict that voids should be filled 
by dwarf galaxies which are not observed (Hoyle & Vogeley 
2004). 

In this paper we are going to present a scenario of void 
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formation which is consistent with astronomical observa¬ 
tions and leads to the formation of low density regions from 
very small initial fluctuations that existed at the moment of 
the last scattering. 


As will be demonstrated, the evolution of voids is a non¬ 
linear process. That is why this approach is based on exact 
solutions of the Einstein equations. 

The structure of this paper is as follows: in Sec. [5] the 
model to be used is presented. In Sec. HU the initial condi¬ 
tions at the last scattering moment are presented. These 
data constrain values of parameters and functions used in 
the theoretical models. Since the influence of the gradient of 
pressure cannot be neglected, in Sec 0] an approximate and 
an analytical solution of a model with inhomogeneous dis¬ 
tribution of radiation is presented. This model qualitatively 
shows that the gradient of pressure plays a significant role 
in the evolution of the Universe after the last scattering mo¬ 
ment. In Sec. |S] results of and fully non linear evolution of 
models with inhomogeneous radiation distribution are pre¬ 
sented. In Sec. IQ constraints on the values of cosmologi¬ 
cal parameters are derived from the fact that the structure 
formation proceeds differently in various cosmological back¬ 
ground models. 
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2 THE SPHERICALLY SYMMETRIC 
INHOMOGENEOUS SPACE-TIME 


A spherically symmetric metric in comoving and syn¬ 
chronous coordinates is of the form: 

ds 2 = e A(t ’ r) c 2 dt 2 -e B< ~ t ’ r) dr 2 -R 2 {t,r) (d6 2 + sin 2 6d(j> 2 ) .( 1 ) 


The Einstein field equations for the spherically sym¬ 
metric perfect fluid distribution (in coordinate components) 
are: 



2 ^~+2A"+ A' 2 -A'B^j=~Kp + k, (5) 


where e is the energy density and p is the pressure. ' stands 
for dt and ' stands for d r . For a mixture of matter and 
radiation, we have: 


Cmat + tradi 

( 6 ) 

1 

^ trad- 

( 7 ) 


The Einstein equations can be reduced to (Lemaitre 
1933): 

nR 2 R!e = 2 M\ ( 8 ) 

KR 2 Rp = -2M, (9) 

where M is defined by: 

2 M(r,t) = R{r,t) + R(r,t)e~ A< ' t ’ r ^ R 2 (r,t) + 

— e~ B( - t ’ r ^R ,2 (r,i)R(r,t) — ^AR 3 (r,t). (10) 

In a Newtonian limit Mc 2 /G is equal to the mass in¬ 
side the shell of radial coordinate r. However, it is not an 
integrated rest mass but active gravitational mass that gen¬ 
erates a gravitational field. As it can be seen in eq. 1101 the 
mass is not constant in time and in the expanding universe 
it decreases, as seen from @. 

From the equations of motion T 01 ^- f 3 = 0 we obtain: 


T 0a - a = 0 
T la ;a = 0 
T 2a ;a = 0 


B + 4 f = 


2 i 

e + p’ 


A' = - 


2 p' 

e +p’ 



( 11 ) 

( 12 ) 

(13) 


T 3 “ Q = 0 =+ ^=0. (14) 

<J(p 

Equations m and m reproduce the well known fact 
that the perfect fluid energy-momentum tensor inherits the 
symmetries of the metric of the space-time. 

The B function can be derived from eq. ©: 


B(r, t ) = B(r,to) + In R! 2 — 
Using eq. ifHt we obtain: 
B(r, t) _ ^' 2 ( r .*) v 


f ,-A'R 

J dt ^' 

*0 


1 + 2 E(r) 


exp / dt 


t 

/• 


2 R(r, t) 


e(r,t) +p(r,t)] R'(r,t) 


~rP{r,t) , 


(15) 


(16) 


where E(r) is an arbitrary function. 


2.1 The Lemaitre—Tolman model 

If the dust equation of state is considered, the above model 
becomes the Lemaitre-Tolman model. As follows from eq. 
iriit and eq. 1161 . if p' = 0 then: 


= 1 , 


R' 


1 + 2 E 

Then the metric 0 becomes: 


ds 2 = e 2 dt 2 — 


1 + 2 E(r) 


dr 2 — R 2 (t, r)dfl 2 , 


(17) 

(18) 


(19) 


where dH 2 = d 6* 2 + sin 2 9d(j> 2 . Because of the signature 
(+, —, —, —), the E(r) function must obey E(r) > — ^. 

The Einstein equations reduce to the following two: 


2 


upc 


2 M' 
R 2 R' ’ 


( 20 ) 



2 E + 


2 M(r) 
R 


+ lArf. 


( 21 ) 


where M(r) is another arbitrary function and k = ^ 2 + 

When R' = 0 and M' ^ 0, the density becomes infi¬ 
nite. This happens at shell crossings. This is an additional 
singularity to the Big Bang that occurs at R = 0, M' ^ 0. 
Shell crossing can be avoided by setting the initial conditions 
appropriately. 

Equation (ED can be solved by simple integration: 


R 

. =c(t-t B (r)), ( 22 ) 

+ i A R 2 

where ts appears as an integration constant, and is an ar¬ 
bitrary function of r. This means that the Big Bang is not 
a single event as in the Friedmann models, but occurs at 
different times at different distances from the origin. 

Thus, the evolution of a Lemaitre-Tolman model is 
determined by three arbitrary functions: E(r), M(r) and 
ts(r). The metric and all the formulae are covariant under 


/ 


dR 


\/ 2E + 2 ~W 
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arbitrary coordinate transformations of the form r = f(r'). 
Using such a transformation, one function can be given a de¬ 
sired form. Therefore the physical initial data for the evolu¬ 
tion of the Lemaitre-Tolman model consist of two arbitrary 
functions. 


2.2 The Friedmann limit 


Assuming homogeneity, the above model becoms the Fried¬ 
mann model. 

In the class of coordinates used here, one can choose the 
radial coordinate as: 


R(r, t) = ra(t), 

Then: 


(23) 

(24) 


M{r,t) = M 0 r 3 , 

2 E(r) = -kr 2 . 
e A(r ’ 4) 1 , 

B(r,t) _^ a 2 

1 — kr 2 ’ 
6 M 


(25) 

(26) 

(27) 

(28) 
(29) 


where a is the scale factor and k is the curvature index of 
the Friedmann models. Then the metric becomes: 

ds 2 = c 2 dt 2 — a(t) 2 — 7-2 S ^ n2 ' (30) 

The condition 3 = 0, reduces to: 


P + 3- (p + — 0- (31) 

a \ cr ) 

Assuming the dust equation of state, eqs. 0 , and m 
become the Friedmann equations: 


2 !=4^ + I Ac2 ’ (32) 

• 2 ,2,1 4 2,1.22 /qo n 

a = —kc + — kc pa + — Ac a (33) 

The Friedmann limit is an essential element in our ap¬ 
proach. As mentioned above, our model of void formation 
describes a single void in an expanding Universe. Far away 
from the origin, the density and velocity distributions tend 
to the values that they would have in a Friedmann model. 
Consequently the values of the time instants (i.e. initial — 
to and final — tf instants) and values of the density and 
velocity fluctuations are calculated with respect to this ho¬ 
mogeneous background. 


3.1 Observational constraints 


The observations of the CMB provide us with the redshift 
of the surface of the last scattering. Once the redshift is 
known, density of matter, density of radiation, temperature 
and pressure can be estimated. The measured value of the 
temperature fluctuations, AT « 10 -5 can be converted into 
density and velocity fluctuations of a baryonic matter at that 
time. The estimated amplitude of the initial fluctuations, 
inside the voids, are (Bolejko, Krasiriski & Hellaby 2005): 

• density fluctuations - A p/p = 5 = —6 ■ 10 -5 

• velocity fluctuations — AV/Vb = v = 2 • 10 ~ 4 


The observed redshift of the CMB is z ~ 1090. At that 
redshift the amplitude of pressure of matter can be esti¬ 
mated from the perfect gas equation of state: 


P 

pc 2 


k B T 

prriHC 2 


nr 10 . 


(34) 


This estimation shows that the pressure of matter after the 
last scattering moment is negligible in the evolution of the 
Universe. 

The contribution of radiation to the energy density can¬ 
not be neglected and the ratio of radiation energy density 
to matter energy density at the last scattering moment is: 


trad _ gT A 
C-mat PC 2 

where a = 


= £i(l + *), (35) 

4(( and a is Stefan-Boltzmann constant; 


aT^ 

o — — cmb an( j Tcmb is the current temperature of 

' Per 

the CMB background; p cr is critical density and is equal 
to (3H 0 )/(8nG). 

At the moment of last scattering, for Qmat = 0.3 this 


is: 


~ 0.2. (36) 

€mat 

As can be seen, this ratio decreases for low redshifts. How¬ 
ever, at the early stage of the evolution after the last scat¬ 
tering moment the radiation should have an influence on 
the evolution of the Universe. This issue will be considered 
futher. 


3.2 Initial fluctuations in the linear approach 

Although the present-day density contrast inside voids is 
large (5 ~ —0.9) and it is unlikely that the linear approxima¬ 
tion might handle it appropriately, the following calculations 
are presented to provide a comparison. Let us assume that 
the Universe is homogeneous with only small perturbations 
of the form: 


3 CONSTRAINTS ON THE INITIAL 
CONDITIONS 

Voids are vast regions in which hardly any galaxies are ob¬ 
served. An average radius of a void is 12 /U 1 Mpc and the 
density contrast inside is —0.94 (Hoyle & Vogeley 2004). 
Such structures must have evolved from small initial fluc¬ 
tuations that existed at the last scattering instant. In this 
section the estimation of the initial conditions is presented. 


a = ab(l + a(r,t)) (37) 

p = p b {l + 5(r,t)). (38) 

Inserting these perturbations into eq. and using eq. 
m we obtain: 

S + 2—6- l K c 2 p b 5 = 0 (39) 

&b 2 

In the simplest case, i.e. the Einstein-de Sitter model, 
this equation reduces to: 
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s 



2 

3c 2 1 2 


S = 


0 , 


and has the analytical solution: 
S = Ci t 2/3 + C 2 t~\ 


(40) 

(41) 


where r is dimensionless time, i.e. r = t/to- 

As can be seen, the first factor describes the growing 
mode, while the second factor describes the decaying mode. 
Assuming that only the growing mode exists, we obtain: 

8 = 5ot 2 ^ 3 , (42) 


where So is the initial density contrast. 

Assuming that the present-day density contrast of voids 
is <S = —0.9, for the Einstein-de Sitter model, t = txj/t^s ~ 
3.5 x 10 4 (where tu is the age of the Universe and tLS is the 
last scattering instant), we obtain: 

5 0 « 8 x 10 -4 . (43) 


However, due to the very large present-day value of the 
density contrast inside voids the linear approximation is in¬ 
adequate. Moreover, the solutions of the formula m have 
neither upper, nor lower, limit. The mathematical structure 
of this equation allows the solution to be of any value, even 
of large negative value. However from the physical point 
of view it is impossible to have a density contrast lower 
than Sum = — 1. Unlike the linear approximation, the real 
density contrast inside the voids should asymptotically ap¬ 
proach — 1. Therefore, this relation is invalid in investigation 
of void formation and an fully non-linear solution must be 
considered. 


3.3 Initial fluctuations in the Lemaitre—Tolman 
model 

Evolution of voids in the Lemaitre-Tolman model was con¬ 
sidered in Bolejko, Krasinski and Hellaby (2005). In that 
paper the evolution of a void from small initial velocity 
and density fluctuations is investigated. The results imply 
that the existence of present-day voids requires the ampli¬ 
tude of the initial density and velocity fluctuations to be of 
S ~ v k- 5 x 10 -3 . 

These results are by one order of magnitude larger than 
the results obtained in the linear approach. Fig. Q presents 
the evolution of the density contrast in the linear approxi¬ 
mation and in the Lemaitre-Tolman model. Both models are 
of the same initial conditions, So = 3 x 10 -3 . As expected, 
the density contrast decreases more slowly in the non- linear 
model than in the linear approach. The difference between 
the results greater than 10% starts at r ~ 500, which cor¬ 
responds to a redshift of z ~ 15. For redshifts larger than 
z = 15 the linear approach reproduces the results with error 
larger than 10%, and for r « 7000 (z m 2.25) the density 
contrast is below —1 which is unphysical. 

Such a disproportion between the theoretical predic¬ 
tion (<5o ~ 10 -3 ) and the observational measurements (So ~ 
10 -5 ) can be interpreted in two ways. The first alternative 
is that the dark matter is made of some weakly interacting 
massive particles, which at the moment of the last scattering 
could have the amplitude of fluctuation of 1CU 3 . The sec¬ 
ond interpretation is that the matter fluctuations at the last 
scattering moment were of 10 -5 and to obtain results consis¬ 
tent with observations some other effects must be taken into 



Mo 


Figure 1. The evolution of the density contrast in the linear 
approach and in the Lemaitre-Tolman model. Both models are 
of the same initial fluctuations, <5o = 3 x 10~ 3 . For r > 7000 the 
linear approach predicts a density contrast lower than —1, which 
is not physical. 


account. As was presented in Sec. im radiation is not neg¬ 
ligible and might influence the evolution of voids. Although 
it seems rather unlikely that introduction of a small amount 
of radiation (temperature fluctuations at the last scattering 
moment were of ~ 10 -5 ) could have a significant impact on 
void evolution, it will be shown below that in reality the radi¬ 
ation component plays an important dynamical role. It will 
be also shown that the linear approximation is inadequate 
here, and that the influence of radiation must be considered 
in the non-linear formalism of relativity. In the next section 
an approximate model with inhomogeneous distribution of 
radiation is presented. 


4 APPROXIMATE SOLUTION 

As long as all fluctuations are very small the linear approach 
should be able to reproduce the proper results. Let us con¬ 
sider an approximate solution in the linear regime. 

In comparison with Sec. on lot us consider the fi ma t = 
1 , Ha = 0 model with a small perturbation of radiation. 
Despite the existence of radiation let us assume that unper¬ 
turbed quantities are the same as in the Einstein-de Sitter 
model: 


k = 0 —> E(r) = 0 


Ro = rah = ra 0 


to) 


2/3 



1 

6nGt 2 


(44) 

(45) 

(46) 

(47) 


Futhermore, as only the pressure fluctuations are of our 
interest let us assume that initial velocity and density fluc¬ 
tuations are equal to zero. Radiation is as follows: 


£ = £<>(1 + 7 ) 


(48) 
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t 1 ~ 2 ft 

6 = 20 POC U 


(49) 


P = ^ = ~fb(l + 7) = b (1 + 7 o exp (~r 2 /L 2 )) ,(50) 

where L 2 is the size of the perturbed region of radius 5 
kpc at the last scattering moment. 70 = 8 • 10~ 5 and is 
the amplitude of radiation fluctuation. Based on eq. 1451 for 
f Imat = 1 the ratio of radiation energy density to matter 
energy density is 1/20 and that is why in eq. 1491 this value 
appears. 

The density contrast is: 


8 = 


(51) 


r(52) 


P~ Po 
Po 

In the Friedmann limit, eq. 1 123) : 

2 6 M 2 6 M 0 f Ro \ 2 kc 2 

KCpo = ^r = 3 ( v ^J = ^ 2 ' 

From eq. nul l and the assumption that there are no 
density and velocity fluctuations: 

In the following, we assume that A = 0 (at the last 
scattering and a few billion years after the A is negligible), 
and A = 0. 

Under these assumptions: 

s=~ 3 ^r~ 


(54) 


6 tt Gt ' 2 


Let us focus on the numerator first. Using eq. llTUl l: 

e -s(t ’ r) f? , 2 (r,t )- 1 (1 + 2 E{r)) 

R 2 V 

exp — / dt 


R 2 


2 R(r,t) 


-p'(r,t) + 


(55) 


e(r,t) +p(r,t)\ R'[r,t) 

\ to 

_J_ _ X_ 

~ R 2 R 2 

Applying the assumption stated at the beginning of this 
section and assuming that: 

Po Po 


poc 2 + (4/3)p 0 Poc 2 ’ 


X = exp I — / dt 


z 

-f 

t O 


2 R(r,t)p'(r,t) 


e(r,t) + p(r,t)\ R'(r,t) 


- 1 


Applying series expansion to the exp: 

1 o 2r , r\ f , z Por~ 8/3 R 
[por- 2 ] R> 


19 r r z C 

^ = 60^°1? exp(— + 2 } J 1 


1 = 


T 

I 9 r r 2 C 

60 2r 7o^ex P (—X^/dr 


por - 8 / 3 2 
[p 0 r- 2 ] 3r 


Close to the origin, r « 0: 


I= JL r 2 l± I drV - 5/3 

180 L 2 


T 

7 d ” 


(56) 

(57) 

(58) 

(59) 

(60) 


T- 1 2 70 /, -2/3\ 

I= 15 r lA {1 - T } 

As can be seen, X has the maximal value, : 

I = — r 2 ltL 

t max — £ 2 

Inserting the above to eq. lIUH : 


5= -3- 


1 

W 


6 tt Gt 2 


o _2 »2 

__ O C t Q 2/3 

20 L 2 a 2 Q T 


(61) 


(62) 


(63) 


To calculate the value of the coefficient, we assume that 
at to = t ls j a 0 = 1. Then: 


5 = -3 x 10 


- 32/3 


(64) 


Comparing this result with eq. nui it is easy to see 
that introducing only pure radiation fluctuation, the results 
are the same as the results obtained in the linear approach 
without any radiation component. The corresponding initial 
density fluctuation in the linear approach would be: 


8 0 = -3 x 10 “ 


(65) 


This proves that radiation is of great importance in 
voids evolution. However to be able to obtain precise predic¬ 
tions, more accurate calculations are needed. The evolution 
of voids within an non linear model is presented in next 
section. 


5 THE ALGORITHM 

The computer algorithm used to calculate the evolution of 
a void was written in Fortran and consisted of the following 
steps. Numerical methods are from Press et al. (1986) and 
Pang (1997). 


5.1 The initial data 

5.1.1 Time instants 

As follows from Sec. |U]to specify the model, one needs to 
know three functions. In this case these function are: den¬ 
sity, velocity and radiation distribution functions, and are 
presented in subsections below. Because it is required that 
at large distace from the origin the model of a void becomes 
an homogeneous Friedman model, these functions are pre¬ 
sented in form of given fluctuations imposed on the homoge¬ 
neous background. All background values are calculated for 
the time instants obtained from the foilwing formula (Pee¬ 
bles 1993): 


t(z) = 


OO 

k! 


dz 


(1 +Z)y/V{z)' 


( 66 ) 


where: 


T>{z) — fl 7 (l + 2 ) + flmat(l + z) + Hk( 1 + z) + Ha , (67) 

Hois the present Hubble constant, 12 a: = 1 — fl 7 — Qmat — TIa- 
The initial instant (to) is set to be at the last scattering 
moment, which took place when z ~ 1089 (Bennett et al. 
2003) and the final instant (tf) when z = 0. 
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5.1.2 The initial density perturbations 

The radial coordinate was defined as the areal radius at the 
moment of last scattering, measured in kiloparsecs: 

r = i = Ro/d = R(r, to)/d 

where d = 1 kpc. 

The initial density fluctuations, imposed on the homo¬ 
geneous background, were defined by functions of the radius 
l, as listed in Tables 0 and [^] and the actual density fluctu¬ 
ations followed from: 


Po(£) = Pbo ■ S(£), 


( 68 ) 


where p b Q is the density of the homogeneous background and 
at the initial instant, and can be expressed as: 

pbO = Slmatpcr(l + z) 3 . (69) 

The mass inside the shell of radius R(r, to) at the initial 
instant, measured in kiloparsecs, was calculated by integrat¬ 
ing eq. 0 : 


M(£,to) — M (0, to) 



(70) 


Since the density distribution has no singularities or zeros 
over extended regions, it was assumed that Imin = 0 and 


M = 0 at i = 0. 


5.1.3 The initial velocity perturbations 

The initial velocity fluctuations, imposed on the homoge¬ 
neous background, were defined by functions of the radius 
£, as listed in Tables0and|2] and the actual velocity fluctu¬ 
ations followed from: 

U(£) = U b o ■ i/(4 (71) 

where U b o is the velocity of the homogeneous background at 
the initial instant, and can be expressed, from eq. I23t , as: 

U b = -R= - to. (72) 

c c 

In the FLRW models the time derivative of the scale 
factor is given by the formula (Peebles 1993): 

a = aH 0 Vv. (73) 

Consequently: 

U = (l + v)-RH 0 Vv. (74) 

c 

In the Lemaitre-Tolman models the proper-time 
derivative of the areal radius R ( U ) is just equal to - R, 
in our case, in consequence of the metric 0: 

U = - e~ A/2 R. (75) 

c 

5.1.4 The radiation perturbations 

In FLRW models the equations of motion reduce to 
(Plebanski & Krasinski 2005): 

dt(ea 3 ) + 3a 2 dp = 0 (76) 


For matter and radiation background (eqs. 0 and 0 ) 
this equation becomes (to simplify the notation, we denote 

Crad — £ and tmat — p). 

(£ + p)a 3 + 3(£ + p)a 2 a + a 2 a£ = 0. (77) 

After the last scattering moment the radiation has not 
been interacting with matter, so we can assume that the 
evolution of matter is the same as in the Universe without 
radiation, i.e.: 



Then from 171 : 


£ = £ »(?) 4 ' 

In the inhomogeneous Universe the radiation can be 
written in the following form: 

£r ad(r,t) = £(t f )(l + z) 4 ((r,t), (78) 

where £ (r, t) is the function which describes the distribution 
and the evolution of radiation. 

According to the current paradigm, after the last 
scattering moment the distribution of radiation has been 
’’frozen”. Consequently the only change in the radiation dis¬ 
tribution is in the amplitude, which decreases with time, as 
the Universe expands: 

C(r,t) = 4>(t)n/(r) (79) 

The crucial problem is to find the form of the 4> function. 
Luckily, the fluctuations of the radiation are very small, of 
amplitude ~ 10~ 5 , therefore, we can assume that the time- 
dependent amplitude of the radiation is the same as in the 
homogeneous background, so: 

4>=l. (80) 

However, this assumption may have to be modified in 
the future if observational data on the distribution of radi¬ 
ation become more detailed. 

Recapitulating: 

trad{£, z(t)) = 4 — Tq mb (1 + z) 4 ■ -/(£)), (81) 

c 

where a is the Stefan-Boltzmann constant, c is the speed 
of light, Tcmb is the current measured temperature of the 
CMB and 7 describes the initial perturbation of radiation, 
as listed in Tables Q and |0 

Solving 1661 with the bisection method the value of 2 
for some instant t can be found. Substituting this value to 
isn the value of energy density of radiation for this instant 
can be calculated. 


5.2 Computing the evolution 

The algorithm of the evolution consisted of the following 
steps: 

(i) From eqs. 171 and 171 the value of R and from 0 
the value of M can be calculated. Then using the predictor- 
corrector method the value of R(r,t + r) and M(r, t + r) in 
the time step of r can be found. We futher denoted them (as 
all the quantities found in this time step) by the subscript 
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(ii) Once R T and M T are known, we can derive t T from 

181 • 

(iii) Then from 0 and ED we derive p T . 

(iv) From e r and p T we can calculate A r by integrating 

eq. 1121 . 

(v) U T can be calculated as follows: 

From eq. iTHjI and i TTTTO we obtain: 




+e 


B(r,t 0 ) Jt 0 dt(A'(r,t)e A < > '’ f )/ 2 !7(r,t)) / R' (r,f) 


(82) 


By solving this equation by the bisection method, for the 
time t = to + r we can calculate U T . 

(vi) Once U T and A T are known, R(r, t) and futher M(r, t) 
can be calculated. 

(vii) We repeat steps 1-6 until t = tf. 


6 RESULTS 


The measurements of density contrast inside voids are based 
on the observations of galaxies inside them (Hoyle & Vogeley 
2004). However, because in central regions no galaxies are 
observed, the real density distribution inside the voids is 
unknown. Assuming that luminous matter is a good tracer 
of dark matter distribution and extrapolating the value of 
the density contrast measured on the edges of voids (where 
galaxies are observed) into the central regions of voids, we 
can conclude that the density inside the voids is below the 
value of 0.06pf,. This value will be called the limiting value. It 
is expected that the model will predict a present-day density 
inside the voids below this limiting value. 

While the problem of choice of the cosmological model 
is open and the efforts to determine the values of f! ma t and 
Ha still continue, in this section (except subsection 16.31 we 
focused only on the cosmological model which is the most 
popular in the cosmological comunity, that is ACDM model. 

Figure |21 shows the shape of the initial perturbations. 
The explicit forms of these perturbations are presented in 
TablesQand|21 The results are presented in FigCU In four out 
of seven model voids were formed. As can be seen, models 
with inhomogeneous distribution of radiation have no prob¬ 
lems with reproducing regions of density below the limiting 
value. 

To compare our results with the observational data, Fig. 
2 ] presents the average density contrast inside the voids as a 
function of a relative distance from the origin. The average 
density contrast was calculated as follows: 

5 = ^4^ - 1, (83) 

P 


where p is the present background density and < p >: 


3 Ale 2 

<p> ~ 4tt GR 3 ’ 


(84) 


Curve 1 presents the results of run 1, as listed in Table 
Q Curves NGP and SGP corespond to density contrasts 
of voids in the 2dFGRS data estimated by Hoyle & Vogeley. 
Although the profiles match at the center they do not fit 
accurately at the edges of voids. In our model the density 
contrast tends to increase faster than the observed one which 


Table 1. The initial Density Perturbations used in the runs. 


Run 

Profile 

Description 

1 

<5_ = 1 — JF 

7 + = l + |^ 
v = 1 + Q 

Isocurvature—like perturbation. 

Reconstructed the present-day voids. 

2 

<5_ = 1 -T 

7— = 1 — 
v = 1 + G 

Adiabatic perturbation. 

The collapse after 20 millions years. 

Leads to high-density regions. 

3 

S- = 1 - T 

70 = 1 
v = 1 + Q 

Isothermal perturbations. 

Do not lead to low-density region. 

4 

5+ = 1 + T 

7+ = 1 + %T 
v = 1 + Q 

Adiabatic perturbation. 

Reconstructed the present-day voids. 

5 

1 + + 

II II H 

1 + + 

<Oco|rfi. V) 

Vi 

Reconstructed the present-day voids, 
although the density fluctuations are 
positive and velocity perturb, are negative. 

6 

<5o = 1 

7+ = 1 + | T 
i/o = l 

Reconstructed the present-day voids. 

7 

<5_ = 1 - T 

7 = 5- = 1 — T 
v = 1 + Q 

Leads to the cluster formation, 
with the value of the central density 

21 • 10 3 p 6 . 


could be caused by too strong assumptions about the evolu¬ 
tion of radiation (O and 111 and would suggest that the 
distribution of radiation did evolve after the last scattering 
moment. There is another possibility that could explain the 
difference between these two profiles. The density contrast 
estmated by Hoyle & Vogeley is based on observations of 
galaxies. It is possible that the existence of dark matter in 
walls would make the real density higher than that estimated 
from counting of galaxies. 


6.1 Initial perturbations 

Introducing radiation into the calculation we need to know 
the relation between matter and radiation perturbations. In 
linear theory there are three concepts of these relations: 

(i) adiabatic perturbations, where 7 = |<5, 

(ii) isocurvature perturbations, where 7 = |(<5 — &>); (80 
is some initial value of 5) 

(iii) isothermal perturbation, where 7 = 0 . 

It should be stressed that in realistic conditions there 
are no pure adiabatic or isocurvature, or isothermal per¬ 
turbations and the relations between density and radiation 
perturbations are more complicated. However, it is instruc¬ 
tive to know what kind of relation is more suitable for the 
process of void formation. 

The results presented in Figure [2] imply that voids can 
be formed out of adiabatic or isocurvature perturbations and 
there is no significant difference between these two forms of 
perturbations, as long as the gradient of the radiation is neg¬ 
ative. With an isothermal perturbation low density regions 
cannot be formed as the gradient of radiation is important 
in the process of void formation. 
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Figure 2. The shape of the initial fluctuations. 



Figure 3. The density distribution inside the void. The numbers 
of curves correspond to the numbers of runs (description in Table 
m in the ACDM background model. The curve denoted by vl 
(void limit) refers to the measured density inside the voids. 



r /Rvoid 

Figure 4. Comparison of the density contrast obtained in run 1 
(1) with the observed density contrast from Hoyle and Vogeley 
(2004) (SGP) and (NGP). 


Table 2. Detailed descriptions of the initial fluctuations from 
Table 0 All the values in the table are dimensionless, and the 
distance parameter is the areal radius in kiloparsecs £ = Ri/ lkpc. 


Profile 

Parameters 

T{t) = a ■ f c -e d f + g 

a = {A — g) ■ {—b)~ c ■ e c 

S = l-b 

g = A ■ [1 — (—b) c ■ e -c ■ ft.]” 1 
h = (R- b) c e d ( b ~ R ) 
c = d ■ b 
d= -0.2 
b = —4.6 

R = 40 

A = 6■10 -5 

Q(£) = A ■ (b • arctan c — d ■ i— 

II 

1 

o 

1 

2 .2 -2 
f ■ e 9 — e h — e 9 — 

b = 4 

m • e~ n ) • k + p 

c = 0.02 • i — 0.02 

d = A 

/ = 0.7 
g = i 

h - f -1 

7 

7= ^ 

J 3 

m = 1.225 

n - £=22 
n 12 
i. — i 

1+0.03^ 

p = -2 x 10 -5 


6.2 Validity of the linear approach 

Fig. 0 presents the evolution of the density contrast in 
three models: in the linear approximation, in the Lemaitre- 
Tolman model, and in an inhomogeneous model with radia¬ 
tion perturbation only. The first two models are of the same 
initial condition. As one can see, the linear approach, given 
by 1391 . is inadequate. At the early stages radiation drives 
the evolution. Because this effect is not considered in the lin¬ 
ear approach, the results obtained in the linear approach are 
not accurate. For later times, when the radiation becomes 
negligible, the density contrast is too large to be correctly 
handled by the linear approximation. 

6.3 Constraints on the background model 

Figure 0 presents the evolution of voids (the initial con¬ 
ditions like in runs 1 and 3) in four different background 
models: 

(a) Qmat = 1, Ov = 0, 0, =4.77- 1CT 5 , 

(b) rimat =0.4, n A = 0, fi 7 =4.77- ltr 5 , 

(c) Qmat =0.27, Ov = 0, a, =4.77- 10“ 5 , 

(d) = 0.27, Ov = 0.73, n 7 = 4.77 • 10" 5 . 

These results imply that in the absence of radiation, or 
of the gradient of radiation, the structure formation goes on 
faster in the models which are filled with a greater amount 
of matter (curves 3a, 3b, 3c and 3d - for more details see 
Bolejko, Krasiriski & Hellaby 2005). Voids cannot be formed 
within this kind of radiation perturbations. By introducing 
a realistic distribution of radiation voids are formed more 
likely. 

According to the astronomical observations, about 40% 
of the volume of the Universe is taken up by voids (Hoyle 
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Figure 5. Comparison of the density contrast evolution in the 
linear regime (upper line), Lemaitre—Tolman model (middle line) 
and the model with only radiation fluctuations. The initial fluc¬ 
tuation in the linear regime and in LT model are <5o = 3 x 10 —3 . 
In the model with only radiation fluctuation, the initial radiation 
perturbations are 70 = 8 x 10 “ 5 . 



Figure 6. The density distribution inside the void with three 
different background models. The numbers of curves corespond 
to the numbers of runs (description in Table ITI and to the back¬ 
ground model (as listed in sec. ioi . The curve denoted by VL 
(void limit) refers to the measured density inside voids. 

and Vogeley 2004). This means that void formation is not 
an isolated event, but it is a very probable process. Thus, 
it can be used to put some constraints on the cosmological 
background model. 

The results in Fig. [3] imply that the presence of radia¬ 
tion is important for void formation. The contribution from 
radiation to the evolution of the system is more significant 
in the models with smaller value of Q m at- Therefore, the 
process of void formation constrains the value of 

As can be seen in Fig. El the cosmological constant is 
not needed to reconstruct the present-day voids, thus the 
process of void formation does not constrain the cosmolog¬ 
ical constant. Although there is a difference between the 
evolution in a model with A and without, these differences 
are small and by choosing different initial conditions can be 


minimised. The model with = 0.4 barely reaches the 

limiting value in the center of the void while models with 
ffmat ~ 0.3 fit the observations best. 


7 CONCLUSIONS 

The aim of this paper was to build a model of void forma¬ 
tion which would simulate the evolution of voids and would 
be consistent with observational constraints. We developed 
a model that describes the process of void formation from 
small initial velocity, density and temperature fluctuations, 
that existed at the moment of the last scattering, and fully 
recovers very low values of density contrast inside the voids. 
However, in our theoretical model, the present density in¬ 
creases faster at the edges of voids than in the observed 
profiles. There could be several explanations: 

(i) Other shapes of the initial perturbations would repro¬ 
duce satisfying results, 

(ii) The assumption (consistent with the widely accepted 
paradigm) that the distribution of radiation did not evolve 
from the last scattering moment is not fulfilled. 

(iii) Matter around voids can distinctly depart from 
spherical symmetry. 

(iv) The real density contrast increases faster than the 
density contrast of luminous matter. 

The main conclusion of this paper is that until sev¬ 
eral million years after the last scattering moment radiation 
cannot be neglected in models of structure formation. The 
gradient of radiation is significant in the process of void 
formation. The negative gradient of radiation causes faster 
expansion of the space inside the void, hence the density 
contrast decreases faster there. The excess of radiation pres¬ 
sure simply drives matter out of the region destined to be a 
void and piles it up on the edges. This effect is purely rela¬ 
tivistic and in the model presented above radiation does not 
interact with matter (this is an accurate assumption for the 
time after the last scattering). As a result to evolve struc¬ 
tures like voids the amplitude of density fluctuation at the 
last scattering moment does not have to be larger than 10 -5 . 
Thus, the fluctuation of dark matter at the last scattering 
can be of the same amplitude as the fluctuation of baryonic 
matter. 

The process of void formation can put some limits on 
the values of cosmological parameters. It was found that 
models with Qmat ~ 0.3 describe the present voids best. 
In models with larger value of the cosmological constant 
the density contrast was lower, due to a longer period of 
evolution. Unfortunately, the only estimation of the density 
distribution inside the voids is based on the observations of 
galaxies, and since in the central parts of voids no galax¬ 
ies are observed, there are no precise estimations of density 
contrast inside them. Therefore we cannot make any precise 
conclusions about the value of the cosmological constant. At 
present the evolution of voids in both CDM models, with and 
without the cosmological constant, is in agreement with the 
observational data. 
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